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1.  WORK  OVERVIEW 


The  first  section  in  this  report  serves  as  a  brief  overview  of  the  work.  The  second  section 
describes  an  analytical  theory  of  the  reflection  and  transmission  of  light  from  nonlin¬ 
ear  optical  interfaces.  The  third  section  presents  an  application  of  the  general  theory  to 
some  simple  nonlinear  waveguide  configurations.  The  last  section  deals  with  a  very  im¬ 
portant  aspect  of  above-threshold  ionization,  namely,  the  intensity  dependence  of  total 
versus  differential  ionization  rates. 

1.1  GENERALIZATION  OF  THE  LAWS  OF  REFRACTION  AND  REFLECTION 
FOR  A  NONLINEAR  INTERFACE 

We  have  generalized  Snell’s  and  Fresnel’s  formulas  of  reflection  and  refraction  for  the 
case  of  a  nonlinear  interface.  We  have  explicitly  dealt  with  the  case  when  a  plane  wave 
is  incident  from  a  linear  medium  onto  an  interface  with  a  medium  of  Kerr-type  nonlin¬ 
earity.  We  have  investigated  an  incident  transverse  electric  (TE)  polarized  wave,  positive 
nonlinearity  constant  in  the  nonlinear  medium,  and  totally  reflecting  interface  at  low  in¬ 
cident  intensity. 

Our  results  can  be  summarized  as  follows.  We  could  show  that  the  solutions  of  the  non¬ 
linear  wave  equation,  describing  stationary  behavior,  fall  into  three  categories: 

1.1.1  Uniform  phase  wave.  The  aimplitude  is  inhomogeneous,  depending  explicitly  on 
the  position  in  the  nonlinear  medium.  This  wave  plays  a  crucial  role  in  the  regime  of 
total  internal  reflection  (TIR). 

1.1.2  Uniform  amplitude  wave.  The  phase  is  that  of  a  plane  wave,  but  with  wave  vector 
depending  on  the  intensity.  This  wave  is  excited  in  the  nonlinear  (NL)  medium  in  the 
entire  regime  outside  that  of  TIR,  and  also  above  a  certain  intensity  threshold  in  the 
regime  of  TIR. 

1.1.3  Inhomogeneous  phase  and  amplitude  waves.  These  waves  play  some  role  in  a  small 
part  of  the  TIR  regime  and  at  extremely  high  incident  intensities.  Therefore,  they  are 
hard  to  observe,  and  nonstationary  phenomena  will  probably  be  highly  competing.  By 
matching  the  above  nonlinear  waves  with  the  linear  incident  and  reflected  waves  at  the 
boundary,  we  managed  to  generalize  Fresnel’s  and  Snell’s  formul^ts  for  the  case  under 
consideration.  Outside  the  regime  of  TIR  (incident  angle  9i  <  9^,  where  9c  is  the  criti¬ 
cal  angle  for  TIR)  the  usual  Fresnel’s  emd  Snell’s  formulas  hold  with  the  index  of  refrac¬ 
tion  in  the  NL  medium  being  the  actual  intensity-dependent  expression.  In  the  regime 
of  TIR,  9c  <  9i  <  7r/2  ,  we  defined  a  second  critical  angle  9o.  In  the  interval  9c  <  9i  < 
9o  TIR  is  maintained  up  to  a  certain,  relatively  low,  incident  intensity  Iq  .  Above  this 
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threshold  TIR  is  lost  and  normal  reflection  and  refraction  takes  place.  The  nonlinea.- 
evanescent  waves  of  1.1.1  go  continuously  over  into  the  NL  radiation  wave  of  1.1.2.  We 
derived  explicit  expressions  for  the  phase  shift  and  reflectivity  associated  with  phenom¬ 
ena  in  this  Oe  <9i  <  9^  interval,  both  for  I  <  Iq  and  I  >  Iq. 

In  the  interval  <  ^''/2,  TIR  is  maintained  up  to  a  relatively  high  incident  inten¬ 

sity  /q.  Above  this  threshold  TIR  is  lost  and,  since  the  inhomogeneous  waves  of  1.1.3  are 
extremely  sensitive  to  small  perturbations,  nonstationary  phenomena  (i.e.,  generation  of 
pulses)  take  place. 

1.2  INTRINSIC  PROPERTIES  OF  NONLINEAR  TE  WAVES  GUIDED  BY  A 
SYMMETRIC  SLAB 

We  have  shown,  by  an  appropriate  scaling  of  parameters  characterizing  the  nonlinear 
TE  guided  modes  of  a  symmetric  slab  (linear  core  embedded  between  two  identical  me¬ 
dia  with  Kerr  nonlinearity),  that  their  mutual  dependence  (dispersion  relations)  can  be 
given  by  particularly  simple  expressions.  This  permits  us  to  determine  the  range  of  vari¬ 
ation  of  these  parameters  and  the  conditions  for  the  appearance  of  asymmetric  waves 
and  surface  modes.  It  also  hints  at  the  possibility  of  bistable  switching  with  hysteresis 
between  different  branches  of  the  dispersion  curves. 

1.3  INTENSITY  DEPENDENCE  OF  TOTAL  VERSUS  DIFFERENTIAL 
IONIZATION  RATES  IN  ABOVE-THRESHOLD  IONIZATION 

Above-threshold  ionization  (ATI)  refers  to  the  ionization  of  an  atom  by  an  intense  laser 
field,  such  that  the  ejected  electron  absorbs  more  laser  photons  than  is  necessary  for  ion¬ 
ization.  The  energy  spectrum  of  the  ejected  electrons  thus  consists  of  separate  peaks 
corresponding  to  the  absorption  of  0,  1,  2, ...  photons  beyond  the  minimum  number. 

The  heights  of  the  peaks,  i.e.,  the  partial  ionization  rates,  are  strongly  dependent  on  the 
intensity  of  the  laser  field.  For  fields  with  an  intensity  >  10‘^  W/cm*  their  intensity  de¬ 
pendence  can  no  longer  be  explained  in  terms  of  perturbation  theory  with  respect  to  the 
laser  field.  In  contrast,  the  total  ionization  rate,  i.e.,  the  sum  of  all  the  partial  rates,  is 
strictly  proportional  to  with  N  the  minimum  number  of  photons  which  have  to  be 
absorbed  for  ionization.  This  striking  behavior  has  been  observed  in  all  pertinent  ex¬ 
periments  up  to  /  >  10**  W/cm^.  It  is  also  exhibited  in  computer  simulations  for  the 
hydrogen  atom.  The  same  phenomenon  has  also  been  observed  (theoretically)  in  nuclear 
beta  decay  in  an  intense  laser  field,  where  the  total  decay  rate  is  unaffected  by  the  field 
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up  to  the  critical  field  strength,  while  the  differential  rates  are  already  dramatically  dis¬ 
torted  by  fields  which  are  less  intense  by  ten  orders  of  magnitude.  In  both  cases,  this 
phenomenon  is  due  to  the  fact  that  the  ionization  or  decay  takes  place  much  like  a  two- 
step  process.  The  first  step  is  the  actual  ionization  or  decay,  while  the  second  step  is  a 
final-state  interaction  between  the  laser  field  and  the  ejected  electron,  in  the  process  of 
which  the  electron  energy  is  redistributed.  We  have  started,  baused  on  previous  experi¬ 
ence  with  nuclear  beta  decay,  to  attempt  to  justify  this  two-step  assumption  from  first 
principles.  In  particular,  we  want  to  find  out  up  to  what  intensity  the  dependence 
will  persist. 

These  specific  areas  of  investigation  are  described  in  detail  in  the  following  sections. 
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2.  GENERALIZATION  OF  THE  LAWS  OF  REFRACTION  AND 
REFLECTION  FOR  A  NONLINEAR  INTERFACE 


2.1  INTRODUCTION 

We  derive  a  wave  equation  describing  stationary  phenomena  along  the  interface.  The 
solution  in  the  linear  medium  is  the  usual  plane-wave  one.  We  show  that  in  the  NL 
medium  the  solutions  fall  into  three  categories:  (a)  waves  with  a  constant  phase  (nonlin¬ 
ear  evanescent  waves),  (b)  waves  with  constant  amplitude  (nonlinear  plane  waves),  amd 
(c)  waves  with  nonuniform  amplitude  and  phase.  From  the  boundary  conditions  (match¬ 
ing  the  tangential  components  of  the  electric  and  magnetic  fields  at  the  interface),  we 
obtain  the  NL  generalization  of  the  Fresnel  formulas.  The  implication  is  the  following. 

Let  Ml  stand  for  the  index  of  refraction  in  the  linear  medium  (ci  =  n’  is  the  dielec¬ 
tric  constant)  and  the  linear  part  of  the  index  of  refraction  in  the  nonlinear  medium 
(cj  =  n]  is  the  linear  dielectric  constant  in  the  NL  medium).  If  >  rii,  then  in  the 
entire  interval  0  <  $  <  ir/2  o(  angles  of  incidence  0  (9  is  measured  from  the  interface 
and  not  from  the  normal  to  it),  the  waves  in  the  NL  medium  are  the  nonlinear  plane 
waves  and  the  phenomenon  is  very  similar  to  the  phenomenon  at  a  linear  interface.  In 
particular,  everywhere  in  Fresnel’s  and  Snell’s  formulas  (j  should  be  replaced  by  the 
intensity-dependent  expression  «2  +  alEjp  (a  nonlinearity  constant,  amplitude  in 
the  NL  medium).  If  n,  <  then  in  the  9^  <  9  <  irf  2  interval  (cos*  9^  =  (7/(1),  the 
situation  is  similar  to  the  1x7  >  ni  case.  In  the  interval  0  <  9  <  9^,  TIR  is  possible  up  to 
a  certain  critical  intensity  /«,  =  Icr{9). 

Our  main  conclusion  is  that  in  this  interval  there  exists  a  second  critical  angle  9o  (sin’  9o 
=  j  sin’  9c  or  approximately  9o  =2  9c /V2),  such  that  in  the  interval  9o  <  9  <  9c  the  non¬ 
linear  evanescent  waves  exist  up  to  a  certain  critical  intensity  Icr ,  and  at  this  intensity 
they  go  continuously  over  into  the  nonlinear  plane  waves;  i.e.,  in  this  interval  the  critical 
intensity  (upper  limit)  of  TIR  coincides  with  the  threshold  intensity  7,^  (lower  limit)  of 
the  nonlinear  plane  waves.  In  the  0  <  9  <  9c  interval,  however,  Icr  >  Ith,  and  in  the 
intensity  range  7,^  <  7  <  7c,  both  waves  are  possible.  Which  one  of  them  will  exist 
at  a  given  intensity  depends  on  the  past  of  the  system  (hysteresis  in  the  reflection  m  a 
function  of  intensity  being  one  possible  consequence).  A  linear  stability  analysis  confirms 
these  findings  and  leaves  no  room  for  the  nonuniform  amplitude  and  phase  waves.  Thus, 
their  existence  requires  further  study. 
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2.2  GEOMETRY,  WAVE  EQUATION,  BOUNDARY  CONDITIONS 


We  deal  with  the  geometry  described  in  the  following. 

A  plane  wave  with  amplitude  E’o  is  incident  on  the  interface  between  a  linear  medium 
(dielectric  constant  fi )  and  a  nonlinear  medium  with  Kerr-t\pe  nonlinearity  (dielectric 
function  Cj  +  a|£l2j*).  The  wave  is  incident  from  the  linear  medium,  its  amplitude  is  fTo 
and  the  angle  of  incidence  is  6.  The  amplitude  of  the  reflected  wave  is  rEo ;  r  is  the  re¬ 
flection  coefficient.  The  incident,  reflected  and  refracted  waves  are  all  polarized  perpen¬ 
dicularly  to  the  plane  of  incidence,  along  the  y-axis  (TE  polarization).  The  amplitude  of 
the  transmitted  wave  is  Ej . 

According  to  the  assumption  about  TE  polarization,  we  look  for  the  electric  field  E  and 
magnetic  field  H  in  the  form 


Here  u  is  the  frequency  and  0  is  the  propagation  constant,  related  to  the  wave  number 
in  vacuum  k  =  a//c,  index  of  refraction  ni  =  y/e[,  and  angle  of  incidence  5  as  = 
k^ei  cos^  0. 

Obviously  div  E  =  0.  From  curl  E  =  —  ^  ^  we  obtain 


H,=-{0lk)E-,  H.  =  -{i/k)-^  (2) 

Using  these  equations  one  can  see  that  div  =  0  is  automatically  satisfied.  From  curl 
H  =  -  ^  we  obtain 

c  at 


-h  kcE  +  0Hi  —  0 


(3) 


with 


€ 


ei ,  if  X  <  0 

ej  -i-a|E2j^  if  X  >  0 


(4) 
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Using  Equation  2  in  Equation  3,  we  can  eliminate  and  H,  and  obtain 


(PE 

dx^ 


[k^€-0^]E  =  0 


which  is  our  starting  equation. 

Boundary  conditions  require  continuity  of  the  tangential  components  of  the  electric  and 
magnetic  fields  at  the  interface,  i  =  0,  i.e.,  continuity  of  E^  and  H, .  In  view  of  Equa¬ 
tions  1  and  2,  this  means 


E  and  E'  continuous  at  i  =  0  (6) 

Equations  5  and  6,  together  with  the  definition  of  k,  c,  and  U  form  the  basis  of  the  the¬ 
ory. 


2.3  SOLUTION  OF  THE  WAVE  EQUATION  IN  THE  DIFFERENT  REGIONS 
2.3.1  Solution  in  the  Linear  Medium 

From  Equation  4  we  find  that  in  the  linear  medium  (i  <  0)e  =  Ci.  Using  this  in 
Equation  5  yields 


E"  -\rk\E  =  0 


(7) 


with 


k\  =  ti  sin^  9  (8) 

ki  is  just  the  x  component  of  the  full  wave  vector  ky/T^.  The  solution  can  be  written  as 
the  sum  of  an  incident  and  a  reflected  plane  wave; 

E(i)  =  E'oe’*'^ -h  rE’oC'**'*,  z<0  (9) 

where  Eo  is  the  amplitude  of  the  incident  plane  wave  and  r  is  the  reflection  coefficient. 
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2.3.2  Solution  in  the  Nonlinear  Medium 


From  Equation  4  we  find  that  in  the  nonlinear  medium  (i  >  0)e  =  £2  —  o(\E\‘ .  Using  this 
in  Equation  5,  we  find 


E"  ■-[k^e^-S'\E-k^-a\E\^E  =  Q  (10) 

In  the  rest  of  this  section  we  categorize  the  possible  solutions  of  this  equation. 

2. 3. 2.1  Uniform  Phase  iVave  (Nonlinear  Evanescent  Wave) 

Assume  that  E  =  £'2(1)6''®  where  £2  is  real.  For  this  case 


£" -y  £2 +«*'£*  =0  (11) 

where  7*  =  /3'*  —  fc’£2-  A  first  integral  can  be  found  by  multiplying  this  equation  with 
£1(=  ^^).  It  is  given  by 


(£;)^  -  ^^El  +  lak^E*  =  const. 


(12) 


For  evanescent  waves  7*  >  0  (i.e.,  costf  >  cos^e  =  y  ^  ^  both  £2  and 

£4  — ►  0  if  X  — ►  00,  hence  const  =  0.  It  is  then  straightforward  to  integrate  Equation  12 
once  again,  yielding 


£2  =  ± 


7  1 

k  cosh[7(i  —  io)| 


or  £  =  £2  e'^ 


(13) 


Here  Xq  is  a  constant  of  integration  to  be  determined  from  the  conditions  of  continuity. 
.\  constant  ph2ise  and  amplitude  wave  (£,  =  E'f  =  0)  can  also  be  found  from  Equation 
11  or  Equation  12: 


aEl  =  €i  cos*  0-62 


(14a) 


or 
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Cl  cos*  $  —  ocEl  +  £2 


(146) 


This  corresponds  to  an  intensity  which  makes  a  given  incident  angle  critical,  the  right- 
hand  side  of  Equation  14b  being  the  square  of  the  effective  index  of  refraction  in  the 
nonlinear  medium. 

Solutions  of  Equation  12  for  an  arbitrary  constant  can  be  expressed  in  the  form  of  ellip¬ 
tic  integrals.  We  come  back  to  this  in  2. 3. 2. 3. 

2. 3. 2. 2.  Uniform  Plane  and  Amplitude  Wave 

Assume  now  that  E  =  £^2^'**’*^'^*  where  E^  =  const.  From  Equation  10  we  find  that 
this  is  possible  if 


fc*  =  fc*  (£2+ a|£;2 !")-/?"  >0  (15) 

If  £2  >  £i,  such  a  wave  always  exists.  If  £2  <  £i,  then  for  £1  cos^  <  £2(d  >  6^)  such 
a  wave  exists  always.  If  £i  cos^  B  >  fj(B  <  Bg),  then  such  a  wave  exists  if  ct\E2\^  > 

£i  cos*  d  ^  £2  =  AhW,  he-)  if  the  intensity  exceeds  a  certain  threshold  intensity.  For 
I  =  Itk  we  just  reobtain  the  constant  phase  and  amplitude  wave  of  Equation  14a.  The 
nonlinear  plane  wave  can  thus  be  written  as 


=  x>0  (16) 

where  is  the  intensity-dependent  wave  vector  as  given  by  Equation  15,  and  £’2  is  to 
be  determined  form  the  conditions  of  continuity. 

2. 3. 2. 3.  Nonuniform  Plane  and  Amplitude  Wave 

Assume  that  both  the  amplitude  and  phase  are  now  functions  of  z,  i.e., 


£(x)  =£2  (17) 

Substituting  this  into  Equation  10  we  can  separate  the  real  and  imaginary  parts  of  the 
wave  equation,  yielding 
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E'^  -  +  {k^e2  -/?"  +  k'^a\E2\^)E2  =  0 


(18) 


from  the  real  part,  and 


2A:j£';  -h  k!,E2  =  0 

from  the  imaginary  part.  Equation  19  can  be  integrated  right  away,  giving 


(19) 


£j  •  fcj  =  const  =  Cl  (20) 

i.e.,  the  expression  on  the  left-hand  side  is  an  integral  of  motion.  The  solution  in  section 
2. 3. 2.1  corresponds  to  =  const  =  0,  and  the  solution  in  section  2. 3. 2.2.  corresponds  to 
both  E2  and  =  const.  Using  Equations  20  and  18  we  can  eliminate  /bj  and  obtain 

£7  +  [k"(7  -0^  ^  k^aEl\E:  =  0  (21) 

^7 

Again,  by  multiplying  Equation  21  with  E^  we  can  find  a  first  integral  in  the  form 

(Kf  +  ^  ;3’|£?  +  \eaZ'.  =  c,  (22) 

Equation  12  is  a  particular  case  when  Ci  =  0  .  The  solution  of  Equation  22  can  again  be 
reduced  to  elliptic  integrals. 

2.4  BOUNDARY  CONDITIONS,  FRESNEL’S  AND  SNELL’S  FORMULAS 
2.4.1  Tht  Case  fj  >  £1 

We  can  match  an  incident  and  reflected  plane  wave  to  a  nonlinear  plane  wave  in  the  en¬ 
tire  0  <  <  7r/2  interval.  Inserting  Equations  9  and  16  into  Equation  6,  we  obtain 


Eo(l  4-  r)  =  £36*'" 


(23) 
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from  the  continuity  of  £  at  i  =  0  and 


k,E,[l-r)=k^E^t'-  (24) 

from  the  continuity  of  E'  at  i  =  0.  By  comparing  the  imaginary  parts  of  Equations  23 
and  24  we  see  that  <p  =  0,  r  =  real.  Using  Equation  23  in  Equation  24,  we  can  eliminate 
E2  from  Equation  24  (note  that  k^  also  depends  on  E2  in  view  of  Equation  15),  yielding 


ki 


1  -r 
1  +r 


—  fcj 


(25) 


If  we  take  the  square  of  this  equation  and  use  Equation  8  for  k^ ,  and  Equation  15  for  k^ 
{E2  is  substituted  from  Equation  23),  then  we  obtain  a  fourth-order  equation  for  r.  It  is, 
however,  linear  in  the  intensity  parameter  I  =  aE^  /ci .  We  can  therefore  express  /  as  a 
function  of  r,  giving 


I 


(r-ro)(r-  ;^)  \ 

(r  +  l)‘  / 


(26) 


with 


ro  =  -cos^e 

sin  fl  -I-  ^/(ca/ei)’  -  cos*  $ 

We  can  now  plot  /  as  a  function  of  r  and  invert  this  plot  graphically  to  obtain 
r  =  r(I,$,e2,ei).  First  we  note  that  if  a  =  0  (linear  media  on  both  sides  of  the  inter¬ 
face)  then  7  =  0,  and  from  Equation  26  r  =  ro  or  r  =  Since  -1  <  Tq  <0,  the  only 
physically  meaningful  solution  is  r  =  Tq  .  This  is  just  Fresnel’s  formula  for  the  reflection 
of  a  TE  wave  from  a  linear  boundary.  Furthermore,  the  r.h.s.  of  Equation  26  is  positive 
if  -1  <  r  <  To  <0,  i.e.,  we  have  to  plot  Equation  26  only  in  this  interval.  The  result  is 
shown  in  Fig.  1.  The  reflection  coefficient  for  low  incident  intensity  coincides  with  that 
given  by  the  (linear)  Fresnel  formula.  With  increasing  intensity  it  decreases  monotoni- 
cally  and  approaches  -1  as  7  -*■  00.  The  interface  becomes  less  and  less  transparent  as 
the  intensity  increases. 
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2.4.2  The  Case  fo  <  Ci 


We  can  define  a  critical  angle  of  incidence  for  which  cos^  6^  =  ^.  In  the  linear  case 
for  0  <  ^  <  total  internal  reflection  takes  place  and  for  f  >  6  >  6c  ordinary  reflection 
takes  place.  In  the  nonlinear  case  we  see  from  Equation  15  that  for  6  >  0^  the  nonlin¬ 
ear  plane  waves  exist  for  any  /,  whereas  for  B  <  9c  they  exist  above  a  certain  threshold 
intensity  Also,  from  Section  2.3,2,  the  nonlinear  evanescent  waves  exist  in  this  inter¬ 
val,  below  a  certain  critical  intensity  4.  Therefore  we  consider  the  9  >  Be  and  9  <  9c 
cases  separately. 

2.4.2.1  9c  <9  <  7r/2 

In  this  case  much  of  the  treatment  of  the  previous  subsection  holds.  In  particular.  Equa¬ 
tion  26  can  be  rewritten  for  this  case  as 


(28) 


where  vq  is  still  given  by  Equation  27.  Now  0  <  <  ro  <  1  and  the  r.h.s.  is 

positive  if  r  <  ro .  Again,  for  7  =  0,  the  solution  r  =  ro  reproduces  Fresnel’s  formula; 
otherwise  we  plot  I  in  the  interval  -1  <  r  <  ro.  The  result  is  shown  in  Fig.  2.  The 
reflection  coefficient  for  low  incident  intensity  coincides  with  that  given  by  the  (linear) 
Fresnel  formula.  With  increasing  intensity  it  decreases  monotonically.  The  intensity  I  = 
sin^  9c  (i.e.,  Cj  -I-  aE^  =  ej  makes  the  interface  transparent,  and  the  reflection  is  0.  If  we 
increase  the  intensity  further,  the  effective  index  of  refraction  becomes  larger  in  the  NL 
medium  than  in  the  lineu  one,  and  the  behavior  is  similar  to  that  found  in  Equation  23 
for  £2  >  £i .  Note  that  there  is  a  phase  change  from  0  to  x  at  this  point. 

2.4.2.2  0  <9  <9c 

In  this  case  the  nonlinear  evanescent  waves  and  nonlinear  plane  waves  may  coexist.  We 
treat  them  again  separately. 

2. 4. 2. 2.1  NL  Plane  Waves  (NLPW) 

If  £2  <  £i  and  9  <  9c,  then  from  Equation  15,  we  see  that  the  NLPW  exists  above  the 
threshold  intensity  7,^  given  by  =  0.  Using  that  E,  =  (l  -I-  r)£’o  and  at  threshold 
r  =  1,  we  obtain  (7,^  =  aE^.^/e^) 


cos^  9  -  cos^  9c 
4 


0  <  ItH 


(29) 
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or 


Cl  cos^  d  =  €2  +  a(2J&o)?h  (30) 

In  other  words,  at  the  threshold  intensity  the  given  incident  angle  tf(0  <  9  <  6^)  becomes 
critical  in  the  sense  that  Ci  cos^  9  matches  the  intensity-dependent  effective  dielectric 
constant  of  the  NL  medium.  If  we  increase  the  intensity  further  while  keeping  9  fixed, 
the  condition  of  total  internal  reflection  cannot  be  satisfied  any  more. 

The  continuity  Equations  23  and  24  remain  the  same  for  this  case  and,  consequently. 
Equation  28  also  holds.  The  only  difference  is  that  ro,  as  given  by  Equation  27,  is  just  a 
phase  factor,  Tq  =  ,  where 


tan^  = 


v/cos*  9  —  COS^  9e 
sin0 


(31) 


Since  these  waves  exist  only  above  a  threshold  intensity,  they  do  not  have  a  linear  limit. 
It  is  therefore  more  convenient  to  write  Equation  28  for  this  case  as 


/  = 


COS^  9  —  cos*  9e 

(1  +  ry 


(1  -  rf 

(l  +  r)« 


sin^  9 


(32) 


In  Fig.  3  we  plot  this  expression  for  a  fixed  ^(0  <  as  a  function  of  r.  The  inter¬ 

pretation  of  /  =  sin*  9e  and  the  phase  change  going  from  I  <  sin*  9c  to  I  >  sin*  9c  is  the 
same  as  in  2.4.2. 1. 


2.4. 2. 2. 2  NL  Evanescent  Waves  (NLEW) 

Using  Equations  9  and  13  in  the  continuity  Equation  6  we  obtain 


^0(1  +  '■) 


from  the  continuity  of  E,  and 


ikiEo{\  -  r) 


>/li 


cosh(7io) 


\  a  k 


tanh(7io) 


cosh(7io ) 


(33) 


(34) 
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from  the  continuity  of  E'  a.t  x  =  0.  These  are  two  complex  equations  to  determine  r 
(which  is,  in  general,  complex),  (p  and  lo  as  a  function  of  S  and  Eq.  Dividing  Equation 
34  by  33  we  obtain 


iki  =  7tanh(T'io) 


(35) 


or,  expressing  r  from  here,  we  get 


r  = 


fct  4- t-)rtanh(7Zo) 
k,,  -  n  tanh(7Zo ) 


(36) 


Hence  r  is  a  phase  factor  with  unit  modulus.  Comparing  the  phases  of  both  sides  of  ei¬ 
ther  Equation  33  or  34,  we  see  that 


r  =  (37) 

That  is,  the  phase  shift  of  the  reflected  wave  is  twice  the  phase  shift  of  the  transmitted 
wave,  and  |r|  =  1,  i.e.,  we  have  TIR  where  these  waves  exist. 

If  we  use  Equation  37  in  Equations  33  and  35,  and  then  eliminate  Xo  from  Equation  33 
with  the  help  of  Equation  35,  then  we  obtain  a  biquadratic  equation  for  tan  <p.  The  solu¬ 
tion  yields 


tan’  <p  — 


sin’  $c  —  2  sin’  9  ±  v^sin*  9^  -  8/ sin’  9 
2  sin’  9 


(38) 


From  this  expression  we  infer  the  existence  of  a  second  critical  angle  9^  defined  as 


sin’  9o  =  -  sin’  9, 

2 

or,  since  usually  9c  is  already  a  small  angle,  to  a  good  approximation 


(39a) 


9o  w  9c/\/2 


{39b) 
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We  deal  with  the  intervals  $o  <  $  <  6^  and  0  <  6  <  9o  separately. 

•  $0  <0  <  Oc 

In  this  case  sin’  -  2  sin’  6  <  0  and  only  the  +  sign  in  front  of  the  square  root  in  Equa¬ 
tion  38  gives  physical  results.  For  /  =  0,  Equation  38  becomes  identical  with  Equation 
31  and  reproduces  the  phase  shift  in  total  reflection  for  the  linear  case.  We  have  two  re¬ 
quirements:  the  expression  under  the  square  root  in  Equation  38  has  to  be  nonnegative, 
which  gives 


I< 


sin*  dc 
8  sin’  9 


=  h 


and  the  entire  Equation  38  has  to  be  nonnegative,  which  gives 


(40) 


/  < 


cos’  9  —  cos^  9, 


(41) 


Since  /j  <  /i,  this  means  that  the  nonlinear  evanescent  waves  in  the  9o  <  9  <9^  interval 
exist  if 


0  <  /  <  (42) 

Furthermore,  since  I2  =  2/,^ ,  where  7,^  is  given  by  Equation  29,  in  the  interval  7,^  < 

7  <  2Ith  the  nonlinear  evanescent  waves  and  nonlinear  plane  waves  coexist.  In  Fig.  4 
we  plot  the  phase  of  r  for  a  fixed  9,  as  a,  function  of  7.  The  reflectivity  |r|  is  up  to  I2. 
When  this  is  viewed  together  with  Fig.  4,  it  indicates  the  possibility  of  a  hysteresis  jump 
as  shown  in  Fig.  7.  One  can  increase  the  intensity  to  Tj  from  0  and  have  TIR.  At  this 
intensity  TIR  is  lost  and  the  reflectivity  drops  by  a  finite  amount  corresponding  to  this 
intensity  in  Fig.  3.  If  we  now  decrease  the  intensity  along  this  line  in  Fig.  3,  then  at 
7  =  the  reflectivity  is  1  and  we  match  the  reflectivity  |r|  =  1.  However,  at  this  point 
(i.e.,  when  TIR  is  restored),  the  phase  jumps  by  a  finite  amount  from  0  to  a  \^lue  given 
by  Fig.  4  at  7  =  7t/, . 

•  0<  9  <  9o 

In  this  case  sin’  9c  >  2  sin’  9  and  both  signs  in  Equation  38  are  possible.  The  solution 
with  the  +  sign  exists  for  0  <  7  <  7i,  with  7i  given  by  Equation  41.  The  solution  with 
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the  -  sign  exists  if  Equation  38  and  its  discriminant  are  nonnegative,  giving  <  I  < 

Since  /j  >  again  the  NLEW’s  and  NLPW’s  coexist  for  <  /  <  A-  In  Fig.  5 
we  plot  the  phase  of  for  a  fixed  0  as  a  function  of  I.  The  subscript  ±  corresponds  to 
the  choice  of  sign  in  Equation  38.  Again,  when  this  figure  is  viewed  together  with  Fig. 

3,  it  indicates  the  possibility  of  a  hysteresis  jump  between  the  +  branch  and  the  NLPW 
when  the  intensity  is  increased.  If  it  is  decreased  along  the  NLPW  branch,  it  matches 
continuously  to  |r+  |,  but  now  the  phase  jumps  at  both  branching  points.  It  does  not 
seem  to  be  easy  to  excite  the  -  branch  of  Fig.  5. 

Figure  6  summarizes  the  existence  conditions  on  the  6  —  1  plane.  The  different  bound¬ 
aries  are  given  by  and  12(6).  In  region  I  the  nonlinear  evanescent  wave 

(NLEW)  of  2.4.2. 2.2  exists.  In  region  II  the  NLEW.,.  (positive  branch)  of  2.4.2.2.2  ex¬ 
ists.  In  region  III  the  NLEW  of  2.4.2. 2. 2  and  the  nonlinear  plane  wave  (NLPW)  of  2. 4. 2. 2. 
coexist.  In  region  FV  the  NLEW+  of  2.4.2.2.2  and  NLPW  of  2.4.2.2.1  coexist.  In  region 
V  the  NLEW+,  NLEW_  (negative  branch)  of  2.4.2.2.2  and  the  NLPW  of  2.4. 2. 2.1  co¬ 
exist.  In  region  VI  only  the  NLPW  of  2.4.2.2.1  exists.  Finally,  in  Fig.  8,  we  plot  the 
critical  angle  as  a  function  of  the  incident  intensity  as  given  by  Equation  30.  In  the 
regions  of  coexistence  (III,  FV,  V  )  one  can  expect  bi-  or  multistability.  Which  of  the 
above  waves  is  actually  stable  requires  further  investigation. 

A  closing  remark:  If,  instead  of  a  coherent  incident  field,  the  interface  is  irradiated  by  an 
incoherent  field,  then  the  incident  and  reflected  fields  do  not  add  in  phase  and  the  inten¬ 
sity  at  total  reflection  is  not  at  the  interface,  but  only  2E^.  This  would  increase  the 
threshold  intensity  for  the  excitation  of  NLPW’s  from  Ith  to  2/,^ ,  which  coincides  with 
the  cutoff  intensity  of  NLEW’s  in  the  6^  <  6  <  9^  region.  Therefore,  our  expectation  is 
that  the  observation  of  a  hysteresis  loop  in  the  intensity-dependent  reflectivity  is  almost 
impossible  for  9o  <  6  <  9^.  On  the  other  hand,  for  0  <  <  0o  the  observation  is  possible, 

since  the  cutoff  intensity  of  the  NLPW+  branch  Ii  >  2/,^. 

2.5  REVIEW 

We  have  defined  a  second  critical  angle  9o  in  the  region  of  total  internal  reflection.  If  the 
angle  of  incidence  9  is  6^  <  9  <  9^,  then  the  evanescent  waves  become  cut  off  at  a  certain 
intensity  and  go  over  into  radiation  waves  in  a  continuous  way.  In  the  0  <  9  <  9q  region 
the  reflectivity  may  exhibit  hysteresis  jumps  as  a  function  of  intensity. 

In  the  next  section,  as  an  application  of  this  general  theory,  we  investigate  nonlinear 
guided  waves  in  a  symmetric  slab  configuration. 
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Fig,  1.  Rsfleclion  coeficieut  r  »et3»s  nonlinear  intenaity  ^f=  “5'','=.®+ al£l»  = 
n*  =  {1,51®  is  the  dielectric  constant  of  the  lurear  medium,  e  j  I  I 

(l.eA  a|EP  «  A. 

rplSVr'^StSdtein^^  plot  corresponds^  the  analytical 
wression  (26)  for  «  =  45- .  /  =  oB?  t.  i»  in  arbitrary  units.) 
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Fig.  2.  Reflection  coefficient  r  versus  nonlinear  intensity  I.  (The  meaning  of  the  pa¬ 
rameters  is  the  same  as  previously,  but  with  Ci  =  (1.6)*, ej  =  (1.5)*.  There 
exists  now  a  critical  angle,  cos*  =  cz/ci  or  0^  as  20*,  such  that  for  6  >  6^ 
normal  reflection  and  for  B  <  total  reflection  takes  place.  The  plot  corre¬ 
sponds  to  the  analytical  expression  (28)  for  B  =  45® .) 
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Fig.  3.  Same  as  Fig.  2.  but  for  $  =  15®,  i.e.,  in  the  total  reflection  interval.  (The 
plot  corresponds  to  the  analytical  expression  (32)  with  5  =  15®.  Note  that 
normal  reflection  occurs  only  above  a  threshold  intensity  1,^  (given  by  (29)) 
This  can  be  said  differently:  the  critical  angle,  as  given  by  (29)  or  (30),  be¬ 
comes  intensity-dependent  and  total  reflection,  which  takes  place  at  low  in¬ 
tensity,  is  replaced  by  ordinary  reflection  and  refraction.) 
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Fig.  4.  Tan^«^  versus  the  nonlinear  intensity  parameter  I.  {4>  is  the  phase  shift  upon 
reflection.  Parameter  values  are  the  same  as  in  Fig.  3.  The  plot  corresponds 
to  Eq.  (38)  with  the  upper  {+)  sign  in  front  of  the  square- root  expression.) 
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Fig.  6.  Phase-plane  of  the  reflection  phenomenon  from  a  nonlinear  interface.  (The 
reflection  can  be  char^terized  by  the  angle  of  incidence  9  and  nonlinear 
intensity  \  is  the  threshold  intensity,  above  which  normal  reflection 
exists.  hjfiUs  the  limit  intensity  below  which  total  reflection  exits.  /,  [6)  is 
the  threshold  intensity;  for  I^[9)  <  I  <  IJ9)  a,  second  branch  of  the  total 
reflection  branch  also  exists  Since  1,^9  <  I, [9]  for  I, ^{9)  <  I  <  1,(9),  the 
total  reflection  coexists  with  the  normal  reflection.) 
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Fig,  7.  Reflection  r  vs.  nonlinear  intensity  I  for  total  reflection  and  normal  reflec¬ 
tion  in  their  region  of  coexistence.  (The  loop  indicates  the  possibility  of 
bistable  behavior.) 
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Fig.  8.  Critical  angle  9  vs.  nonlinear  intensity  I.  (The  plot  corresponds  to  Eq.  (30).) 
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3.  ON  THE  INTRINSIC  PROPERTIES  OF  NONLINEAR 
TE  WAVES  GUIDED  BY  A  SYMMETRIC  SLAB 


3.1  INTRODUCTION 

In  recent  years  considerable  theoretical  effort  has  been  devoted  to  the  investigation  of 
nonlinear  optical  waveguides.  This  interest  is  motivated  by  their  capability  to  support 
soliton-type  wave  propagation  and  also  to  exhibit  sudden  switches  between  regimes  with 
different  guiding  mechanisms.  These  features  have  a  number  of  potential  applications. 
Optical  systems,  in  general,  provide  an  accessible  tool  to  study  phase-transition-like  phe¬ 
nomena  occurring  in  nonequilibrium  and/or  nonlinear  systems.  A  typical  device  is  a 
resonator  filled  with  a  nonlinear  medium  (Ref.  1).  In  this  case  the  transmission  of  the 
system  becomes  intensity  dependent,  because  of  the  dependence  of  the  effective  optical 
length  and  absorption  on  intensity. 

The  purpose  of  this  report  is  to  analyze  the  optical  wavequide  bounded  by  nonlinear 
media  as  an  alternative  optical  system  having  many  of  the  described  features.  In  pre¬ 
vious  works,  attention  was  mainly  focused  on  the  problem  of  finding  appropriate  exter¬ 
nal  (control)  parameters  and  on  the  behavior  of  the  system  as  these  control  parameters 
are  varied.  Although  this  is  an  important  issue,  from  an  experimental  point  of  view,  in¬ 
trinsic  features  also  have  to  be  studied  to  achieve  a  full  understanding  of  their  behavior. 
Therefore,  particular  emphasis  is  placed  here  on  revealing  the  mutual  dependence  of  in¬ 
ternal  parameters  which  characterize  a  nonlinear  guided  wave  propagating  in  the  system. 
This  enables  us  to  determine  the  range  of  variation  of  the  parameters,  as  well  as  to  es¬ 
tablish  existence  criteria  for  the  various  waves. 

Historically,  hysteresis  reflection  and  refraction  from  a  single  nonlinear  interface  was  first 
studied  by  Kaplan  (Ref.  2).  The  evanescent  wave  solution  given  there  wcis  later  success¬ 
fully  applied  to  the  study  of  nonlinear  surface  polaritons  (Refs.  3-4).  After  the  initial 
proposals  on  optical  bistability  in  nonlinear  waveguides  (Refs.  5-6),  systems  with  two 
or  more  nonlinear  interfaces  were  studied  in  detail  on  the  basis  of  coupled  mode  theory 
(Ref.  7),  which  is  valid  only  for  weak  nonlinearity.  Application  of  exact  methods  to  the 
nonlinear  slab  followed  thereafter  (Refs.  8-12).  In  an  interesting  paper  Akhmediev  pre¬ 
dicted  bistable  switching  between  different  branches  of  a  given  TE  mode.  Based  on  that 
work  we  elaborate  the  dependence  on  the  propagation  constant  of  various  parameters 
appearing  in  the  solution.  We  show,  in  particular,  that  this  mutual  dependence  of  var¬ 
ious  internal  parameters  suggests  how  these  nonlinear  waves  can  be  excited,  and  deter¬ 
mines  the  associated  field  patterns  completely. 
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3.2  PARAMETER  DEPENDENCE  AND  ASSOCIATED  FIELD  PATTERNS 


Figure  9  shows  the  geometry  of  the  problem  under  study.  A  slab  with  linear  dielec¬ 
tric  constant  c  =  6i  is  embedded  between  two  identical  media  with  Kerr  nonlinearity, 
e  =  -t-  QjE(x),^.  For  TE  polarization  is  the  only  nonvanishing  component  of  the 

electric  field.  It  is  common  to  cissume  the  form  Ey  =  E(x)  exp[i(/?z  -  wt)],  where  /?  is  the 
propagation  constant. 

Then  the  evanescent  wave  solution  for  guided  modes  in  regions  2  and  3  can  be  written  as 
(Ref.  8) 


E(z)  =  ±(2/a)>  (7/A:)sech[7(i-ij.3)];  li|  >  d  (43) 

Here  we  assume  ot  >  0.  We  introduce  the  notation  7’  =  and  ik  =  w/c.  7  is 

uniquely  related  to  0.  It  is  a  constant  associated  with  the  propagation  (attenuation)  of 
evanescent  waves  along  x.  Xj  and  Z3  are  integration  constants  in  regions  2  and  3,  respec¬ 
tively. 

The  possible  wave  patterns  in  the  linear  region  1  are  given  by 


r  cos[«;i(z  -  do)l;  A 


E(z)  =  A 


sin[Ki(z  -  do)]:  5 

- 

cosh[«:3(z  -  do)|;  C;  |z|  <  d 


V  sinh(K2(z  -  do)]:  D 


(44) 


Here  kI  =  fc^£i  —  >0  and  /c*  =  Ki(«3)  is  the  constant  associated 

with  the  propagation  (attenuation)  of  guided  waves  along  z  in  region  1  for  <  0  < 
<  0).  A  and  do  are  further  constants  of  integration. 

Boundary  conditions  require  that  £'(z)  and  E'{x)  be  continuous  at  z  =  ±d.  This  yields 
four  relationships  among  the  parameters  in  Equations  43  and  44.  On  the  other  hand, 
these  field  expressions  contain  five  parameters:  A.Zj.Zjdo,  and  0.  This  was  not  a  prob¬ 
lem  in  the  linear  case:  due  to  linearity,  A  cancelled  from  the  continuity  equations  and 
the  remaining  four  parameters  were  determined  completely  (Ref.  13).  In  the  linear 
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waveguide,  in  particular,  do  —  0,  i.e.,  all  the  field  patterns  are  either  completely  sym¬ 
metric  or  completely  2mtlsymmetric.  Furthermore,  I3  =  -X3  =  —00,  so  that  evanescent 
waves  are  purely  exponentially  decaying  ones.  Also,  /?  is  restricted  to  certain  discrete 
values  (roots  of  the  eigenvalue  equations)  in  the  interval  <  0  <  ky/T[,  so  that  no 

surface  TE  modes  (hyperbolic  patterns  in  1)  exist  for  0  >  ky/e^. 

In  the  nonlinear  case  (a  ^  0)  the  situation  changes  drastically.  0  becomes  an  interval 
continuous  parameter;  i.e.,  instead  of  assuming  discrete  values  in  the  ky/t^  <  0  <  ky/T^ 
interval,  as  in  the  linear  case,  0  can  change  continuously  within  discrete  subintervais  of 
the  above  interval.  Also,  values  in  the  region  0  >  ky/^  (surface  modes)  become  permis¬ 
sible.  This  can  be  seen  from  the  following  treatment.  With  the  help  of  the  four  continu¬ 
ity  equations  we  can  express  any  of  the  other  four  parameters  as  a  function  of  0.  Each 
of  these  equations  can,  with  equal  rights,  be  regarded  as  a  generalized  eigenvalue  equa¬ 
tion  of  the  nonlinear  slab.  They  determine  the  physical  domain  of  parameters,  as  well  as 
the  existence  conditions  for  the  different  types  of  nonlinear  waves  (odd,  even,  asymmet¬ 
ric  guided  modes  and  surface  modes). 

For  later  convenience,  and  in  order  to  bring  the  eigenvalue  equations  to  their  simplest 
form,  it  is  useful  at  this  point  to  introduce  the  scaled  quantities  ij=  x^/d,  x^=  x^jd^ 
do  =  do  /d.  These  dimensionless  parameters  determine  the  extremum  positions  of  the 
field  intensity  distribution  inside  and  outside  the  slab.  We  also  introduce  a?  — 
which  is  a  nonlinear  intensity  parameter  proportional  to  the  intensity  maximum  in  the 
slab.  Instead  of  the  propagation  constant  0  we  define  two  related  parameters:  b]  = 

{k^€i  -  0^)<P  for  <  0  <  ky/(^  and  6^  =  (/?*  -  k0ii)(P  for  0  >  ky/t^.  Finally, 
we  introduce  the  parameter  w’  =  («!  -  €3)Jk’d’,  which  characterizes  the  waveguide  it¬ 
self.  Once  the  waveguide  is  specified,  v  is  fixed.  The  advantage  of  these  scaled  variables 
is  that  seemingly  very  different  waveguides  having  the  same  parameter,  u,  exhibit  similar 
behavior  in  terms  of  the  scaled  parameters.  Thus,  the  usage  of  these  parameters  reveals 
universal  properties  of  the  nonlinear  slab  waveguide. 

We  now  proceed  to  the  investigation  of  the  eigenvalue  equations  and  the  field  patterns. 
We  investigate  the  different  possibilities  enumerated  in  Equation  44  separately.  Further¬ 
more,  we  analyze  in  detail  only  the  cases  when  the  cos  or  cosh  solutions  hold  in  the  slab. 
A  'omplete  treatment  of  all  possible  cases  is  left  to  a  separate  publication  (Ref.  14). 

3.2.1  We  begin  our  investigation  with  the  case  when  the  first  line  holds  in  Equation  44. 
Then  the  boundary  conditions  allow  for  do  =  0  and  do  7^  0,  as  well. 

3. 2. 1.1  If  do  =  0,  the  field  pattern  is  symmetric.  In  this  case,  from  the  boundary  condi¬ 
tions  we  obtain 
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I3=  -  Xj=  1  -  (v^  -  6^)"^''*  tanh“^[6i  (v^  -  taniij  (45) 

which  is  our  first  eigenvalue  equation. 

Instead  of  a*  we  use  the  related  parameter  S  =  cos*  6i ,  which  is  the  nonlinear  inten¬ 
sity  at  the  boundary.  The  other  eigenvalue  equation  then  reads  as 


S  =  u* 


cos*  6i 


(46) 


Equations  45  and  46,  together  with  do  =  0,  characterize  the  symmetric  cosine  wave  com¬ 
pletely,  in  the  sense  that  the  remaining  parameters  are  expressed  as  functions  of  (i.e., 
the  propagation  constant  /?).  Where  in  the  interval  <  0  <  <  &i  <  v) 

this  wave  may  exist  can  be  determined  from  the  following  argument.  The  criteria  of  ex¬ 
istence  Me  simply  that  S  >  0  and  that  the  argument  of  tanh~  ‘  in  Equation  45  should  be 
between  -1  and  1.  Both  these  requirements  lead  to 


b]  <  w*  cos*  bi  (47) 

i.e.,  such  waves  exist  in  those  parts  of  the  0  <  <  v  interval  where  Equation  47 

holds.  As  expected,  all  solutions  of  Equation  47  lie  in  this  interval.  For  any  value  of  v 
there  is  always  at  least  one  such  subinterval  of  bi,  i.e.,  the  lowest  order  symmetric  co¬ 
sine  mode  never  becomes  cut  off.  As  v  increases,  additional  subintervals  may  appear.  In 
these  subintervals  Xs  changes  from  — oo  to  oo  as  is  changed  from  one  bound2U7  of  the 
subinterval  to  the  other.  5  is  zero  at  the  boundary  of  the  subinterval,  reaches  a  maxi¬ 
mum  inside  and  then  falls  back  to  zero  again  at  the  other  boundary.  It  is  interesting  to 
note  that  at  the  right  boundary  of  each  subinterval  the  argument  of  tanh"  ‘  in  Equation 
45  is  equal  to  1,  so  that 

(u*  -  6*)~ tan6i  =  1  (48) 

This  is  just  the  eigenvalue  equation  for  even  TE  modes  of  a  linear  waveguide.  Also,  at 
these  points,  5  =  0.  This  suggests  the  following  mechanism  for  the  excitation  of  these 
nonlinear  waves.  At  very  small  intensity  S  the  system  is  essentially  linear,  and  a  guided 
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wave  which  is  similar  to  a  linear  eigenmode  can  be  launched.  If  we  slowly  increase  S 
while  allowing  bi  (/?)  to  change  according  to  Equation  46,  then  Zj  increases  from  — oo  to- 
w2Lrds  +00.  At  a  threshold  value  of  5  =  Si{Si  =  or  bi  =  nn  for  the  n-th 

subinterval)  we  obtain  X2=  1,  and  a  real  maximum  of  the  field  distribution  outside  the 
slab  appears.  At  5  =  So,»„  >  5i  we  reach  the  maximum  of  Equation  46,  and  5  cannot 
be  increased  further  along  this  branch  of  the  dispersion  curve.  Thus,  this  mode  becomes 
cut  off  at  5  =  5n,,,  (or  at  ii.ma*  corresponding  to  5n,,x»  where  6i,m*x  =  -cotan 
The  other  side  of  this  branch  6i  <  6i  ,o,,x  cannot  be  reached  with  continuous  excita¬ 
tion  from  the  linear  regime.  It  corresponds  to  two  evanescent  waves,  excited  initially 
infinitely  far  away  from  the  slab  boundaries.  Increasing  the  intensity  5  at  the  bound¬ 
ary  attracts  these  waves  closer  to  the  boundaries.  Hence,  they  can  only  be  excited  under 
conditions  very  different  from  those  of  a  waveguide.  Figure  10  shows  the  dispersion  re¬ 
lation  (plot  of  Equation  46)  for  v  =  In  this  case  TEo  and  TE^  modes  are  allowed. 
Also  shown  are  two  typical  TE^  field  patterns  for  S  <  Si  and  Si  <  S  <  Smmx-  The 
guiding  mechanism  for  5  >  5^  is  obviously  different  from  the  usual  total  internal  reflec¬ 
tion  (TIR)  of  linear  waveguides.  It  can  adequately  be  described  as  self-focusing  in  the 
nonlinear  medium. 

This  is  not  the  whole  story,  however.  The  effective  refractive  indexes  on  the  two  sides  of 
the  slab  boundary  become  equal  if  the  intensity  satisfies  €i  =  +  a|£7(d)P  or,  in  our 

notation,  S  =  y-  intensity  the  boundary  becomes  transparent,  and  we  might 

expect  something  strange  to  happen.  That  this  is  indeed  the  case  is  made  clear  in 
the  next  subsection. 

3.2.1.2  Here  we  investigate  the  dg#  0  case,  where  the  field  pattern  is  asymmetric.  The 
analytic  expression  for  dg,  resulting  from  the  continuity  equations,  is 


in^bi  do  = 


1/2 


(49) 


Here  r  =  cos*  bi .  In  this  case  —  I3 ,  and  these  quamtities  are  given  by 


X2.3=  ±1  T  (u*  —  tanh"^  |6i(t;*  —  6J)~ tan  6i(lTdg)  | 


(50) 
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The  nonlineu  intensities  at  the  boundaries  ±d,  S  =±=  a*(6i(l±  do)],  can  be  written  in 
the  more  symmetric  form 


_ 

^  ^  I cos^[6i (1- do)]  cos*[6i(H- do)] 

The  criterion  of  existence  of  these  asymmetric  guided  waves  is  that  0  <  sin*  6i  do<  1, 
or  the  argument  of  tanh"^  in  Equation  50  should  be  between  -1  and  1,  or  S±  >  0.  All 
these  requirements  lead  to 


6?  <  y  cos^  (52) 

i.e.,  such  waves  exist  in  those  parts  of  the  0  <  &,  <  vy/2  <  0  <  in¬ 

terval  where  Equation  52  holds.  For  any  value  of  v  there  is  always  at  least  one  such 
subinterval  of  6i,  i.e.,  the  lowest  order  asymmetric  mode  never  becomes  cut  off.  As  v 
increases,  additional  subintervals  may  appear.  In  these  subintervals,  do  chwges  such 
that  0  <  6i  do  <  x/2,  and  it  is  0  at  both  boundaries  of  the  subinterval.  The  bound¬ 
aries  of  these  subintervals  are  defined  by  the  equality  sign  in  Equation  52.  The  nonlin¬ 
ear  intensity  at  these  points  is  5+  =  Sj  =  S,h  =  v*/2,  as  can  be  seen  from  Equation 
51.  Thus,  the  asymmetric  waves  branch  away  from  the  symmetric  ones  at  the  intensity 
which  makes  the  boundary  transparent.  Consequently,  their  appearance  is  a  threshold 
phenomenon. 

In  Fig.  11  we  have  plotted  the  dispersion  relation.  Equation  51,  and  two  typical  asym¬ 
metric  field  patterns  developing  from  the  symmetric  TEt  mode  for  S  >  S,h.  In  that 
part  of  the  hi  interval  where  the  symmetric  and  asymmetric  waves  coexist,  the  asyz\- 
metric  one  carries  less  energy.  This  may  indicate  that  the  asymmetric  wave  is  generally 
more  stable  and  that  the  part  of  the  dispersion  curve  of  the  symmetric  modes  above  S(h 
is  difficult  to  reach. 

This  dispersion  curve  together  with  the  dispersion  curve  of  the  symmetric  TE2  mode 
of  Fig.  10  hints  at  the  following  interesting  possibility.  If  the  intensity  at  the  boundary 
is  increased  along  S+ ,  then  at  Sm.i  we  reach  the  maximum  of  the  5+  branch.  Denote 
the  corresponding  hi  value  by  .  If  we  try  to  increase  the  intensity  further,  the  sys¬ 
tem  may  jump  over  to  the  5_  branch,  along  which  the  intensity  can  be  increased  up  to 
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the  value  Sn,2-  At  this  value  this  type  of  asymmetric  mode  becomes  cut  off.  If  we  lower 
the  intensity  along  the  5_  branch  we  reach  the  minimum  at  a  point  where  bi  =  6*/^ . 
This  is  the  same  value  where  the  S+  — *■  S_  jump  occurs,  but  with  the  role  of  the  upper 
and  lower  boundaries  interchanged.  The  system  at  this  point  may  jump  over  to  the  sym¬ 
metric  TE2  branch  of  Fig.  10.  (Actually,  these  two  figures  should  be  viewed  together. 
Only  for  clarity  are  they  represented  separately.)  If  now  the  intensity  is  increased  along 
this  branch,  we  eventually  reach  5,i,.  Further  increase  of  the  intensity  excites  the  S+ 
branch  of  the  dispersion  curve.  The  excitation  loop  completed  this  way  would  then  con¬ 
tain  bistable  jumps  with  hysteresis.  Since  the  aisymmetric  waves  carry  less  power  than 
the  symmetric  ones,  lowering  the  intensity  along  5_  requires  less  input  power  than  in¬ 
creasing  it  along  T .  This  is  in  good  qualitative  agreement  with  the  experimental  ob¬ 
servation  of  Ref.  15.  It  is  also  clear  that  this  type  of  bistability  cannot  occur  in  the  low¬ 
est  order  mode,  since  5+  and  S.  branching  away  from  TEo  (the  left  portion  of  Fig.  11) 
do  not  exhibit  this  structure.  This  conclusion  is  also  supported  by  the  experiment. 

3.2.2  Except  that  the  v  dependence  plays  a  more  important  role,  the  antisymmetric  sine 
wave  obtained  when  the  second  line  holds  in  Equation  44  has  similar  properties  to  the 
ones  discussed  above.  Therefore,  we  do  not  elaborate  further  on  this  case.  A  complete 
discussion  is  given  in  Ref.  14. 


3.2.3  We  next  investigate  the  case  when  the  third  line  holds  in  Equation  44,  the  case  of 
the  cosh  wave.  Since  the  boundary  conditions  now  allow  do  =  0  only,  the  wave  is  com¬ 
pletely  symmetric.  The  analytic  expressions  for  the  other  parameters  au’e 


13=  -  Xj  -hi  +  (62  +  v’)  ^^^tanh  ‘  62  (63 -h  v’)"  tanh&2 


(53) 


For  the  nonlinear  intensity  at  the  boundary,  S  =  cosh^  62  >  have 


5  =  u^ 


-h 


cosh’  6* 


(54) 


The  criterion  of  existence  is  S  >  ,  or  the  argument  of  tanh"  ^  in  Equation  53  should  be 

less  than  1,  which  is  fulfilled  in  the  entire  b,  >  O{0  >  k^/TH)  interval  for  all  v.  The  inten¬ 
sity  for  62  =  0  (/?  =  ky/eTj  is  u’,  and  this  value  coincides  with  the  intensity  of  the  lowest 
order  symmetric  cosine  mode  at  61  =  0  (/?  =  ky/T^.  This  suggests  that  the  symmet¬ 
ric  surface  modes  described  by  the  hyperbolic  cosh  solution  can  be  reached  by  increasing 
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the  intensity  of  the  lowest  order  symmetric  cosine  mode  and  tuning  from  (3  <  ky/r[  to 
3  >  ky/r[.  The  appearance  of  the  surface  mode  is  a  threshold  phenomenon  requiring  the 
minimum  nonlinear  intensity  5th  =  i^igure  12  shows  a  plot  of  the  dispersion  relation, 
Equation  54,  and  a  typical  field  pattern.  If  62  is  increased,  5  increases  first,  reaches  a 
maximum  5  =  S^a.^  at  a  certain  62  =  i2.maxt  and  then  decreases  monotonically.  Obvi¬ 
ously,  only  that  part  of  the  dispersion  curve  for  which  S  <  can  be  reached  by  con¬ 
tinuous  excitation  from  the  /?  <  regime.  At  this  point  this  surface  mode  becomes 

cut  off.  The  decreasing  part  corresponds  to  a  very  different  excitation  condition:  surface 
modes  are  excited  outside  the  slab.  It  is  interesting  to  observe  the  behavior  of  12  as  a 
function  of  62-  As  62  is  increased,  r,  increases  first,  starting  from  the  value  1  at  6^  =  0, 
reaches  a  mJiximum  at  62  =  ij.max  and  decreases  again  towards  12=  2  as  62  — ►  00.  That 
is,  the  position  of  the  maximum  of  the  field  distribution  always  remains  in  the  vicinity  of 
the  slab  boundary,  justifying  the  name  “surface  mode.” 

3.2.4  The  case  when  the  fourth  line  holds  in  Equation  44  (hyperbolic  sine  wave)  is  dis¬ 
cussed  in  Ref.  14. 

3.3  CONCLUSION 

The  symmetric  slab  waveguide  bounded  by  nonlinear  media  can  support  symmetric,  an¬ 
tisymmetric  and  asymmetric  guided  waves  and  surface  modes.  The  latter  two  have  no 
counterpart  in  the  linear  case.  Their  appearance  is  a  threshold  phenomenon:  a  thresh¬ 
old  intensity  at  the  slab  boundary  is  needed  for  their  excitation.  By  appropriate  scaling 
of  the  parameters  characterizing  the  different  waves,  we  cast  the  dispersion  relations  of 
these  waves  into  particularly  simple  form.  This  enables  us  to  establish  existence  condi¬ 
tions  for  various  types  of  waves  and  to  determine  the  range  of  variation  of  all  the  pa¬ 
rameters  involved.  Also,  the  interpretation  of  the  results  is  facilitated.  For  example,  the 
threshold  intensity  for  the  appearance  of  asymmetric  waves  is  the  intensity  which  makes 
the  boundary  transparent.  The  threshold  intensity  for  the  appearance  of  a  surface  mode 
is  the  intensity  where  the  lowest  order  guided  mode  terminates.  The  dispersion  relations 
in  this  scaled  form  are  also  very  suggestive  as  to  the  excitation  of  these  waves.  At  very 
small  values  of  the  nonlinear  intensity  at  the  boundary,  the  dispersion  relations  reduce 
to  the  eigenvalue  equations  of  a  linear  waveguide.  It  is  known  how  an  eigenmode  of  a 
linear  waveguide  can  be  launched.  Upon  increasing  the  nonlinear  intensity  adiabatically, 
one  can  consecutively  excite  diflFerent  portions  of  that  branch  of  the  dispersion  relation 
which  originates  from  a  given  linear  limit.  Furthermore,  our  results  seem  to  indicate  the 
possibility  of  a  bistable  loop  in  the  excitation  process  of  higher  order  modes  in  qualita¬ 
tive  accord  with  recent  experimental  findings. 
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Fit;.  1^-  l^ispcrsioti  curve  (intensity  at  the  slab  boundary  vs.  i)roi)agation  constant)  for  tin-  syiniiu  t- 
ric  cosine  modes  ami  for  v  --  5.  (Also  shown  arc  field  |)atterns  belonging  to  the  .S’,  •  an 
‘‘’h. »,  >  ^  |)ortion  of  the  TK  branch  ol  the  <lisi>ersion  curve.) 


Fig.  11.  Dispersion  curve  and  field  pattern  of  the  asymmetric  cosine  wave. 


.  Dispersion  curve  and  6eld  pattern  of  the  symmetric  hyperbolic  cosine  surface  mode 


4.  INTENSITY  DEPENDENCE  OF  TOTAL  VERSUS  DIFFERENTIAL 
IONIZATION  RATES  IN  ABOVE-THRESHOLD  IONIZATION  (REF.  16) 


In  above-threshold  ionization,  total  ion  yields  conform  with  lowest  order  perturbation 
theory  up  to  much  higher  intensities  than  the  individual  peaks  of  the  electron  spectrum. 
This  fact  is  explained  and  an  estimate  of  the  intensity  is  derived  for  which  the  total  ion¬ 
ization  rate  ceases  to  follow  lowest  order  perturbation  theory.  The  derivation  uses  argu¬ 
ments  based  on  simple  dimensional  analysis. 

Multiphoton  ionization  of  atoms  by  intense  laser  fields  has  been  investigated  for  many 
years  (Ref.  1).  In  all  cases,  for  intensities  up  to  about  10^®IV/cm*,  the  total  ionization 
rates  below  saturation  as  measured  via  the  total  ion  yields  were  found  to  be  propor¬ 
tional  to  the  Nth  power  of  the  laser  intensity  I,  with  N  the  minimum  number  of  laser 
photons  which  the  atom  has  to  absorb  to  be  ionized.  This  -  behavior  is  readily  ex¬ 
plained  on  the  basis  of  lowest  order  perturbation  theory  with  respect  to  the  laser  field. 
On  the  other  hand,  the  more  recent  experiments  recording  the  energy  spectra  of  the 
ejected  electrons  start  revealing  markedly  nonperturbative  features  for  intensities  as  low 
as  10^®  IV/cm’  (Refs.  17-22).  A  typical  ATI  (above-threshold  ionization)  electron  spec¬ 
trum  exhibits  equally  spaced  peaks,  separated  by  the  energy  ho*  of  a  laser  photon,  at 
energies 

—  Co  +  nhw  (n  =  0,1,2...)  (55) 

While  perturbation  theory  suggests  that  the  height  of  the  peak  be  proportional  to 
7^  experiments  show  that  for  intensities  a  *m  10*®W/cm*  this  is  no  longer  so:  the 
lowest  peak  (n  =  0)  is  no  longer  the  dominant  one  (peak  switching);  moreover,  the  low- 
order  peaks  become,  for  increasing  intensities,  one  at  a  time  completely  suppressed  (peak 
suppression).  The  heights  of  the  individual  peaks  approximately  scale  with  7^~  >  (Ref. 
17)  for  most  n  (the  lowest  peak  with  a  power  less  than  this)  rather  then  with  7^'*'’*. 
Deng  and  Eberly  (Ref.  23)  identified  an  intensity  I,at  ~  10*^W/cm’  such  that  for  7 
^  7«a«  lowest  order  perturbation  theory  ceases  to  be  valid  for  the  individual  peaks.  At 
the  same  time,  it  still  holds  perfectly  well  for  the  total  ionization  rates  (Ref.  24).  The 
same  behavior  has  also  been  observed  in  numerical  analyses  (Ref.  25).  It  is  the  pur¬ 
pose  in  this  report  to  make  this  apparent  puzzle  plausible  and  to  identify  an  intensity 
hat  »  hat  so  that  for  7  ~  hat  lowest  order  perturbation  theory  (LOPT)  becomes  in¬ 
applicable  for  the  total  rate.  We  will  proceed  on  the  b2usis  of  simple  dimensional  analysis 
without  assuming  any  specific  model. 
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We  first  notice  that  the  behavior  just  described  appears  to  be  a  quite  common  phe¬ 
nomenon  in  multiphoton  physics.  We  may  summarize  it  by  stating  that  often  the  total 
transition  rate  for  some  process  in  the  presence  of  a  laser  field  can  be  calculated  by  per¬ 
turbation  theory  for  much  higher  intensities  than  the  corresponding  differential  transi¬ 
tion  rates.  For  example,  for  the  case  of  the  free-electron  laser,  the  probability  that 
the  electron  emits  n  laser  photons  while  it  travels  through  the  undulator  is  a  very  com¬ 
plicated  function  of  the  laser  intensity  and  only  calculable  by  LOPT  for  extremely  low 
intensities.  On  the  other  hand,  the  gain,  which  is  proportional  to  np„,  can  be 

obtained  from  LOPT  up  to  intensities  which  are  higher  by  about  eight  orders  of  magni¬ 
tude  (Ref.  26).  Another  well-known  example  is  the  Kroll- Watson  sum  rule  for  electron- 
potential  scattering  in  the  presence  of  a  laser  field,  which  says  that  under  certain  condi¬ 
tions  the  total  cross  section  a  =  (^<»  being  tbe  cross  section  with  accompa¬ 

nying  emission  of  n  laser  photons)  Just  equals  the  cross  section  Oj,,,  in  the  absence  of 
the  laser  field  (Ref.  27).  In  this  case  LOPT  obviously  yields  a  —  Ojr,,  while  it  is  not 
applicable  for  the  - 


Before  returning  to  ATI  we  shall  extensively  consider  laser  enhancement  of  nuclear  0- 
decay  as  an  additional  example.  This  is  a  particularly  convenient  case  because  the  decay 
rate  has  been  evaluated,  and  LOPT  has  been  proven  to  hold.  The  very  simple  result  for 
the  field-induced  enhancement  of  the  decay  is  (Ref.  28) 


P  = 


-To 


(eJg)»h» 

mej 


(56) 


for  E  <  Ec  and  Co  <■  .  Here  E  is  the  field  strength  of  the  (circularly  polarized) 

laser  field,  Ec  =  ~  1.6  x  10‘*V/cm  the  so-called  critical  field,  m  the  elec¬ 

tron  mass  and  Cq  the  nuclear  energy  release  of  the  decay  minux  the  electron’s  rest  mass 
energy.  The  quantity  k  is  a  numerical  factor  (k  =  35/64  for  allowed  /9-decay,  and  k  = 
105/32  for  a  first-order  forbidden  decay).  In  contrast  to  the  simple  form  (Eq.  56)  of  the 
enhancement  of  the  total  rate,  the  differential  decay  rate  with  respect  to  the  energy  E^ 
of  the  emitted  electron,  viz.,  dT ficnldEgy  has  a  very  complicated  dependence  on  the 
field  strength  (Ref.  29)  and  cannot  be  obtained  from  perturbation  theory  except  when 
E/Ec  is  extremely  small. 

We  will  now  try  to  rederive  Equation  56  from  dimensional  considerations.  Let  us  first 
list  the  relevant  variables:  besides  m,e,h,  and  e,  there  is  the  field  strength  E  and  fre¬ 
quency  <jj  of  the  laser  field,  the  nuclear  energy  release  Co  and  the  extent  R  of  the  nuclear 
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wave  function.  A  complete  set  of  dimensionless  quantities  that  can  be  formed  out  of 
these  is  given  by  E/E„eclmc^ yy/me^Rlh,a  =  /he,  and  eE/mue.  We  now  make  the 
following  assumptions  or  observatioxis: 

(a)  A  glance  at  the  form  of  the  matrix  element  that  has  to  be  calculated  to  obtain 
the  exact  result  (Eq.  56)  makes  clear  (Ref.  28)  that  the  enhancement  will  de¬ 
pend  on  even  powers  of  eE  and  that  this  is  the  only  occurrence  of  the  charge  e. 
There  is  then  no  additional  dependence  on  a  =  e’/he; 

(b)  the  enhancement  is  a  quanttim  effect,  so  that  we  should  have  p  — »  0  for  h  0; 

(c)  since  normally  Co/mc^  <<  1,  and  hu/me*  «  1,  the  leading  contribution  to  p 
is  nonrelativistic,  i.e.,  independent  of  c; 

(d)  the  dynamics  are  quasistatic,  i.e.,  the  enhancement  depends  only  on  the  field 
strength  E  and  not  explicitly  on  the  frequency  u  of  the  laser  field. 

By  assumption  (d),  any  dependence  of  p  on  the  vector  potential  a  ==  eE/w  or  on  the  di¬ 
mensionless  parameter  eE/tnue  is  ruled  out.  This  leaves  us  with  the  three  dimensionless 
quantities  E/E„eo/me^  ,y/nu^R/h  to  build  the  dimensionless  enhancement  p.  Obser¬ 
vation  (a)  forces  the  leading  contribution  to  p  to  be  proportional  to  {E/EcY ,  (c)  then 
makes  necessary  the  factor  of  (me^  /£o )’  in  order  to  cancel  the  dependence  on  e.  This 
yields  the  result  (Eq.  56)  which  satisfies  assumption  (b).  The  latter  would  still  be  sat¬ 
isfied  with  an  additional  factor  of  y/meoR/h;  however,  y/rm^R/h  ^  0(10“*),  so  this 
would  not  be  the  leading  contribution.  Assumption  (b)  is  also  satisfied  when  the  r.h.8. 
of  Equation  56  is  multiplied  by  arbitrary  inverse  powers  of  y/me^R/h.  However,  the  en¬ 
hancement  p  would  then  become  proportional  to  inverse  powers  of  the  nucleu  radius  R, 
This  is  highly  implausible  in  view  of  the  fact  that  the  matter-field  coupling  can  be  taken 
as  r  •  E.  We  thus  see  that  the  assumptions  (a)  to  (d)  essentially  prescribe  the  form  (Eq. 
56)  for  the  leading  contribution  to  the  enhancement.  Moreover,  assumption  (a)  is  really 
an  observation.  The  same  applies  to  assumption  (b),  although  in  a  less  obvious  way:  if 
one  describes  the  /3-decay  by  the  quartic  Fermi-interaction,  then  the  laser  field  couples 
to  the  charged  particles  taking  part  in  the  process,  notably  the  electron.  This  has  two 
effects:  first,  a  classical  acceleration  of  the  electron  after  the  decay  has  taken  place.  This 
distorts  the  differential  decay  rate,  but  leaves  the  total  rate  unaffected;  second,  the  decay 
rate  itself  is  altered.  This  comes  from  off-shell  contributions  of  the  electron  propagator 
adjacent  to  the  four-Fermi-vertex.  Hence,  any  enhancement  of  the  decay  must  be  quain- 
tum  mechanical  and  proportional  to  some  positive  power  of  h.  Assumption  (c)  is  justi¬ 
fied  by  noticing  that  relativistic  corrections  would  have  to  involve  positive  powers  of  the 
quantity  fo /me*,  which  is  small  for  the  small  nuclear  energy  releases  we  are  considering. 
(For  Co /me*  ~  1,  assumption  (c)  and  Equation  56  do  not  hold.) 
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The  quasistatic  assumption  (d)  is  crucial  for  the  above  analysis.  It  is  highly  plausible 
for  nuclear  ;9-decay,  which  proceeds  on  a  time  scale  much  faster  than  the  time  scale  u}~  '■ 
of  the  laser  field  and  on  a  spatial  scale  much  smaller  than  the  laser  wavelength.  The  di¬ 
mensionless  quantity  ea/me^  =  tEjmuc  (the  occurrence  of  which  in  the  total  decay  rate 
is  ruled  out  by  assumption  (d))  plays,  however,  the  dominant  role  in  the  interaction  be¬ 
tween  the  electron  and  the  laser  field  after  the  decay  and,  consequently,  in  the  differen¬ 
tial  decay  rate.  This  is  obvious  by  noticing  that  the  solution  of  the  nonrelativistic  equa¬ 
tion  of  motion  mx  —  eEainut  is  v/c  =  ~(eE/muc)coB(jt  -1-  const.  Hence  the  parameter 
cE/muc  governs  the  large-scale  motion  of  the  electron  in  the  laser  field  (which  is  of  the 
order  of  the  wavelength). 

We  now  apply  the  same  reasoning  to  above-threshold  ionization.  For  low  enough  intensi¬ 
ties  the  ionization  rate  per  unit  time  allows  for  a  power  series  expansion  in  terms  of  the 
intensity  I 

r  =  2r,  =ro(i-h2r„/ro)  (57) 

f»aO  ns  1 

where  r„  with  N  as  defined  above.  The  question  then,  which  we  shall  try  to 

answer,  is  for  which  intensity  the  leading  correction  Fj  /Fq  to  the  /^-behavior  will 
become  important.  (For  the  individual  peaks  this  intensity  is  given  by  /«««•)  The  prob¬ 
lem  is  almost  analogous  to  field-enhanced  0-decay.  The  most  notable  difference  is  that 
^-decay  takes  place  also  in  the  absence  of  the  field,  whereas  in  ATI  Fq  ~  .  If  we  de¬ 

scribe  the  atom  by  just  a  radius  R  and  an  ionization  potential  V  (replacing  Cq)  the  list 
of  parameters  is  the  same  as  above,  although  some  of  their  values  are  different.  Conse¬ 
quently,  the  set  of  dimensionless  quantities  is  the  same.  We  again  make  the  assumptions 
(a)  to  (d).  Assumption  (a),  provable  by  inspection  in  the  case  of  /?-decay,  will  hold  for 
ATI  as  long  as  no  additional  energy  or  length  scales  enter  the  problem,  so  that  all  occur¬ 
rences  of  the  elementary  charge  which  are  of  atomic  origin  are  hidden  in  the  constants 
R  and  V .  This  will  not  be  so  if  there  are  near-resonant  intermediate  states.  Even  then, 
however,  these  may  well  predominantly  affect  Fq  and  much  less  so  the  ratio  Fi/Fq.  As¬ 
sumption  (b)  appears  to  be  debatable  since,  after  all,  ionization  can  follow  completely 
classical,  chaotic  routes  (see,  e.g..  Ref.  30).  However,  we  will  see  below  that  this  hardly 
affects  our  estimate  (as  long  as  we  do  not  expand  our  list  of  relevant  parameters).  There 
is  no  problem  with  assumption  (c).  Assumption  (d)  is  less  well  satisfied  for  ATI  than  for 
/7-decay,  but  should  still  be  safe,  since  the  atomic  length  scale  is  much  smaller  than  the 
laser  wavelength. 
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It  is  then  obvious  that  we  obtain  the  same  result  as  before  for  the  leading  correction  to 
the  total  ionization  rate,  viz., 


To  mV^ 


(58) 


with  a  different  undetermined  constant  k'.  Assuming  0(A:')  =  1,  we  have  0(ri/ro)  =  1 
for 


/  \3  _  -Tiat 

^  Ib 


(59) 


where  Eb  =  e/a^{ao  =  h^/me^)  is  the  Bohr  field  and  Ib  =  eEg  /4jr  =  7.1  x  j ew? 

the  corresponding  intensity.  For  V  ~  lOeV  this  yields  7,,,  5  x  10“®/b  ^  j cm? . 

This  defines  the  intensity  7..,  >>  mentioned  in  the  Introduction.  Equation  59  is 
the  maun  result  of  this  section.  It  explains  the  observed  persistence  of  the  I"  behavior 
of  the  total  ionization  rate,  as  opposed  to  the  nonperturbative  behavior  of  the  ionization 
rate  into  individual  peaks,  and  predicts  an  order  of  magnitude  for  the  intensity  for  which 
corrections  will  become  importamt. 

Assumption  (b)  could  be  relaoced  to  allow  for  Fi/Fo  to  be  independent  of  K  or  even  pro¬ 
portional  to  negative  powers  of  A  as  long  as  Ft  is  proportional  to  a  positive  power.  The 
r.h.s.  of  Equation  58  woiild  then  be  multiplied  with  some  power  of  the  dimensionless 
quamtity  y/mVR/h.  However,  for  V  =  lOeV  and  R  =  10" 'em  we  have  y/mVRfh  =  1.15, 
so  that  the  order  of  magnitude  of  Ft  /Fq  is  unchanged.  If  Ft  /To  is  considered  to  be  re¬ 
lated  to  an  (Af  +  1)  -  quantum  process,  the  expression  given  in  Equation  58  seems  to  be 
the  more  logical  one.  It  can  readily  be  checked  that  the  leading  relativistic  corrections  to 
Ft/Fo  are  very  small. 

It  should  be  emphasized  again  that  our  assumption  (d)  implies  that  the  total  ionization 
rate  is  independent  of  the  quantity  eE/mue  =  ea/mc*.  This  latter  quantity  character¬ 
izes  the  ponderomotive  potential,  which  has  a  very  significant  impact  on  the  ATI  elec¬ 
tron  spectrum  for  intensities  I  ~  /.«,  (Ref.  31).  Since,  however,  the  former  is  related 
to  the  large-scale  quivering  motion  of  an  electron  in  a  laser  field,  it  should  play  no  role 
in  the  total  ionization  rate.  Hence,  even  when  many  of  the  low  energy  pe2dcs  of  the  elec¬ 
tron  spectrum  are  already  suppressed,  the  total  rate  still  is  proportional  to  .  For  in¬ 
tensities  up  to  lO^^W/em*,  this  has  been  experimentally  observed  in  some  cases  (Ref. 
24).  As  has  been  remarked  by  Lompr4  et  al.  (Ref.  22),  it  may  therefore  be  misleading 
to  think  of  the  ponderomotive  potential  as  added  to  the  ionization  potential. 
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Of  course,  the  predictive  power  of  a  result  derived  by  dimensional  arguments  is  limited. 
Should  the  experimental  value  for  7„,  turn  out  to  be  in  major  disagreement  with  Equa¬ 
tion  (59),  this  could  be  due  to  (1)  a  value  of  the  numerical  factor  it'  significantly  differ¬ 
ent  from  unity,  (2)  more  parameters  than  assumed  above  playing  a  decisive  role  or  (3) 
any  of  the  assumptions  (a)  to  (d)  being  violated.  Any  of  these  findings  would  be  of  in¬ 
terest  by  themselves.  In  any  event,  since  an  experimental  value  is  presently  not  known, 
Equation  59  gives  an  estimate  of  where  one  might  expect  deviations  from  the  be¬ 
havior. 
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